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An ab-initio description of atomic nuclei that solves the nuclear many-body prob- 
lem for realistic nuclear forces is expected to possess a high degree of predictive 
power. In this contribution we treat the main obstacle, namely the short-ranged 
repulsive and tensor correlations induced by the realistic nucleon-nucleon inter- 
action, by means of a unitary correlation operator. This correlator applied to 
uncorrelated many-body states imprints short-ranged correlations that cannot be 
described by product states. When applied to an observable it induces the cor- 
relations into the operator, creating for example a correlated Hamiltonian suited 
for Slater determinants. Adding to the correlated realistic interaction a correction 
for three-body effects, consisting of a momentum-dependent central and spin-orbit 
two-body potential we obtain an effective interaction that is successfully used for 
all nuclei up to mass 60. Various results are shown. 



1. Introduction 

In the last years exact ah initio calculations of light nuclei have become 
feasible with Greens Function Monte Carlo calculations^ and in the No- 
Core Shell Model ^. Here realistic interactions that fit the nucleon-nucleon 
scattering data and the deuteron properties are used Additional three- 
body forces are needed and are adjusted to the spectra of nuclei. Chiral 
perturbation promises to provide a consistent derivation of two- and three- 
body forces 

1 
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2. The Unitary Correlation Operator Method (UCOM) 

Our aim is to perform ab initio calculations of larger nuclei with realistic 
interactions like the Bonn or Argonne potentials in a Hartree-Fock picture 
or a many-body approach with configuration mixing. 

The repulsive core and the strong tensor force of the nuclear interac- 
tion induce strong short-range radial and tensor correlations in the nu- 
clear many-body system. These correlations are in the relative coordinates 

= f'i — f'j and thus can not be represented by products of single-particle 
states like Slater determinants 

\^)=A{\qi)^...®\qA)] (1) 

that are usually used as many-body states in Hartree-Fock or a shell-model 
calculations. A denotes the antisymmetrization operator and | qi ) the 
single-particle states. 

Instead we treat the radial and tensor correlations explicitly by a unitary 
correlation operator C that acts on uncorrelated product states | ^ ) 

I ^ > = C I ^ ) (2) 

such that the many-body state | ^ ) contains the short ranged correlations. 
For the correlator we make the following ansatz 

C = C^'C^ =exp|-i^^f^^jj -expl-i^^^^jj . (3) 

i<j i<j 

It is the product of a radial correlator and a tensor correlator C^, both, 
expressed with a hermitian two-body generator in the exponent. 

2.1. Cluster expansion 

As the ansatz for the correlator contains a two-body operator in the expo- 
nent any correlated operator will contain many-body parts. For example a 
Hamiltonian consisting of one- and two-body parts will turn into 

H = C^HC = C^{^Ti + Y, Vij)C (4) 

= E^^+E^^?+ E ^;S + --+E^i'+ E ySi+-- > 

i i<j i<j<k i<j i<j<k 

where the upper script ["^^ indicates irreducible n-body operators. Here we 
introduce an approximation by keeping terms only up to two-body opera- 
tors. This approximation should be good for systems where the range of 
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the correlator {grij = and g^ij = for r^j > Rc) is short compared to 
the mean particle distances. In that case the probability to find 3 particles 
simultaneously within the correlation range Rc is small. 

2.2. Radial correlator 

The radial correlator (described in detail in shifts a pair of particles 
in the radial direction away from each other so that they get out of the 
range of the repulsive core. To perform the radial shifts the generator of 
the radial correlator uses the radial momentum operator pr together with 
a shift function s(r) that depends on the distance of the two nucleons. 

9rij = - {Prijs{rij) + s{rij)prij^ (5) 

The shift function s{r) is optimized to the potential under consideration. 
It is large for short distances and will vanish at large distances. 

The effect of the transformation | ^ ) ^ | ^ ) is shown in the upper 

(2) 

part of Fig. 1 where the two-body density ^ is displayed as a function of 

the distance vector (ri— between two nucleons in ^He. On the l.h.s. p^pj. 
is calculated with the shell- model state | ) that is just a product of 

4 Gaussians. It has a maximum at zero distance which is in contradiction to 
the short ranged repulsion of the interaction. This inconsistency is removed 
by application of the radial correlator Cr that moves density out of the 
region where the potential is repulsive. The corresponding kinetic, potential 
and total energies are displayed in the lower part of the figure for three 
nuclei. The radially correlated kinetic energy ( C},TCr ) increases somewhat 
compared to ( T ) but this is overcompensated by the gain of about -25 MeV 
per particle in the correlated potential energy. Nevertheless the nuclei are 
still unbound. 

2.3. Tensor correlator 

The tensor force in the S=l channels of the nuclear interaction depends on 
the spins and the spatial orientation f = (ri— r2)/(|ri— r2|) of the nucleons 
via the tensor operator 

S,2{r, r) = 3(0-1 • f)(o-2 • f) - (0-1 • 0-2) = 2 (3{S • rf - S") . (6) 

An alignment of f with the direction of the total spin S = |(cri + cr2) 
is favored energetically. The tensor correlator Cq = exp< —'^^i^jgnij r, 
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Figure 1. Upper part: Two-body density p^glp{ri — r2) of ^He for a pair of nucleons 
with isospin T = and parallel spins, S = Ms = 1. Bright ares denote large probabilities. 
Arrow indicates spin direction and (x,y,z) = (ri — r2) relative distance vector. Lower 
part: corresponding kinetic, potential and total energies per particle of ^He, -"^^O and 
^^Ca, without, with radial, and with radial and tensor correlations (Bonn- A potential). 



defined by the generator 

gnij = ^{rij)^(^{(TiPnij){(Tjrij) + {(Tirij){ajp^^j)^ , (7) 

achieves this alignment by shifts perpendicular to the relative orientation 
Tij. For that the generator of the tensor correlator uses a tensor operator 
constructed with the orbital part of the relative momentum operator = 
Pij—Pr ij ' The r-dependent strength and the range of the tensor correlations 
is controlled by i^{r). For details see ^. 

The application of the tensor correlator Cq leads to the two-body den- 
sity depicted in the right hand contour plot of Fig. 1. One may visualize 
the action of Cq as a displacement of probability density from the 'equator' 
to both 'poles', where the spin of the 5' = ! component of the nucleon pair 
defines the 'south-north' direction. Again this costs kinetic energy but now 
the many-body state is in accord with the tensor interaction and one gains 
the binding needed to end up with about -8 MeV per particle. 
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bare potential 



correlated interaction 

Figure 2. Bare and correlated Argonne V8' interaction in momentum-space. Matrix 
elements range from -2 to +2 fm"-*^ on the l.h.s. and from -1.5 to fm"-*^ on the r.h.s. 

3. Interaction in momentum- space 

The inclusion of the short-range correlations achieves a pre-diagonalization 
of the nuclear hamiltonian that is illustrated in Fig. 2. On the l.h.s. the 
effect of the radial correlations is shown in the channel. The correlated 
interaction evaluated in momentum-space is more attractive and does not 
possess the large off-diagonal matrix elements of the bare interaction. Also 
the tensor components of the correlated Hamiltonian do not connect to high 
momenta as is illustrated with the matrix elements between the and 
the ^1^1 channel. The correlated interaction is therefore a low- momentum 
interaction very similar to the Viow-k ^• 

4. Effective interaction 

To test the two-body approximation we performed no-core shell model cal- 
culations with the correlated AV8' interaction for ^He and compared with 
exact results i^'^i. It turned out that neglecting the n = 3- and 4-body 
parts T^"^] and V^'^'^ of the correlated Hamiltonian leads to an overbind- 
ing which is of the same order as the contribution from genuine 3-body 
forces. Because of the low- momentum nature of the correlated Hamilto- 
nian =T+T[^] + F[^] it can be used directly with simple model spaces 
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Figure 3. Deviation of mean- field binding energies from measured ones. 

built of Slater determinants. The effects of missing higher-order contribu- 
tions to the correlated interaction and of genuine three-body forces is for 
now effectively described by a momentum- dependent central and spin-orbit 
correction Hcorr with four parameters adjusted to four doubly magic nuclei. 
The resulting effective interaction H^ff =H^^ + Hcorr is used for all nuclei 
up to mass number 60. The expectation value of Hcorr is typically 15% of 
the correlated interaction energy. 

5. Hartree-Fock calculations 

For the nuclei listed in the nuclear chart Fig. 3 we minimized the expectation 
value ( ^ I Heff — Tcm | ^ ) / ( ^ | ^ ) with respect to all parameters of the 
single-particle states 

<x|g) = ^c,exp|-^^^||Xi)®|0 (8) 

that are contained in the Fermionic Molecular Dynamics state | ^ ) 
(Eq. (1)). The summation is either for one or for two independent Gaus- 
sians per single-particle state. 

The inclusion of a second Gaussian improves the masses of p-shell nuclei 
substantially (inset of Fig. 3). The largest deviations occur for nuclei with 
an a-cluster structure like ^Be or ^^C or for intrinsically deformed nuclei at 
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Figure 4. One-body densities of intrinsic states. Upper part: PAV calculation -"^^Be, 
20Ne, 28Si (local minimum); lower part: VAP calculation ^^c, i^q, ^^Qs.. 

the middle of the sd-sheh. Some examples of intrinsic shapes are displayed 
in the upper half of Fig. 4. In ^"^Be the extra six neutrons have pulled 
together the well localized pair of a's that form ^Be. The peanut like shape 
for ^^Si is a local minimum. 

6. Projection after variation and variation after projection 

To improve the many-body Hilbert space we project on spin and parity after 
variation (PAV). We also perform variation after projection (VAP) calcula- 
tions in the sense of the generator coordinate method. The intrinsic state 
is minimized here with constraints on radius, dipole moment, quadrupole 
or octupole moment. 

With the above described effective interaction ^^O gains in a VAP cal- 
culation about 5 MeV in binding by forming a-clusters (see lower part of 
Fig. 4) compared to the spherical closed-shell configuration which repre- 
sents the energy minimum in the PAV case. Even for ^^Ca VAP leads to 
a deformation towards a tetrahedron of ten a's which however are much 
more amalgamated than in ^^O (Fig. 4). 

In a next step the different intrinsic shapes obtained in the VAP process 
can be used to perform multiconfiguration calculations. 
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Figure 5. One-body densities of intrinsic ^'^C states. 

A very interesting nucleus is ^^C for which the shell- model configuration 
(l<5i/2)^ (1^3/2)^ is competing with the 3-a cluster structure i^'^^. Like 
in -'^^O the PAV ground state of ^^C turns out to be a spherical shell- 
model state (Fig. 5, V/PAV). However, variation after projection leads to a 
triangular shape made of three a's (Fig. 5, VAP) and 7.2 MeV of additional 
binding as indicated in Table 1. A pure a-cluster configuration obtained 
in a VAP calculation has larger distances between the a-clusters and is 
4.3 MeV less bound than the shell-model configuration. The difference to 
the full VAP is due to the polarization of the a-clusters. The description 
of this polarization is significantly improved by using two Gaussians per 
single-particle state. 

The description of the excited states requires an enlarged Hilbert space. 
Multiconfiguration calculations with the configurations shown in the lower 
part of Fig. 5 show only a small increase in binding for the groundstate 
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Figure 6. Calculated and experimental level scheme for ^"^C. 



but have a significant effect on the radius and the B{E2) value for the 
0^ — > 2^ transition. The additional configurations have been chosen to 
give lowest energies for the 3f and the second and third 0+ states. We find 
3 a-cluster structures for the O2" state and ^Be + a-cluster structures for 
the O3" state. For an improved description of these states a larger number 
of configurations with greater distances between the a's is needed. This is 
consistent with an assumed Bose condensed state^^ for the O2" state. 
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